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We show how to construct a multi-qubit control gate on a quantum register of an arbitrary 
size A'^. This gate performs a single-qubit operation on a specific qubit conditioned by the 
state of other A'^ — 1 qubits. We provide an algorithm how to build up an array of networks 
consisting of single-qubit rotations and multi-qubit control-NOT gates for the synthesis of 
an arbitrary entangled quantum state of A*' qubits. We illustrate the algorithm on a system 
of cold trapped ions. This example illuminates the efficiency of the direct implementation of 
' the multi-qubit CNOT gate compared to its decomposition into a network of two-qubit CNOT gates. 
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I. INTRODUCTION 



■ Entanglement is probably the most intriguing aspect of quantum theory [Q. It attracts due attention not only for 
its epistemological importance Q but also as an essential resource for quantum information processing. In particular, 
quantum computation |^], quantum teleportation quantum dense coding certain types of quantum key 
^ ' distributions ||8|] and quantum secret sharing protocols^, are rooted in the existence of quantum entanglement. 
— , Recently, lot of progress has been achieved in investigation of various properties and possible application of quantum 
\^ ' entanglement. Nevertheless, many questions are still opened. In particular, it is the problem of multi-particle 
I entanglement [ p^ . Specifically, in contrast to classical correlations, quantum entanglement cannot be freely shared 
. among many objects ||ri|, It has been shown recently |lj] that in a finite system of qubits with N{N — l)/2 
entangled pairs the maximal possible concurrence (a specific measure of entanglement |l^, |l^) is equal to 2/N. This 
, value of the bipartite concurrence is achieved when the TV qubits are prepared in a totally symmetric state |S), such 
that all except one qubit are in the state |1), i.e. 



|0), = ^ (lOll ... 1) + 1101 ... 1) + 1110 . . . 1) + • • • + 1111 . . . 0) j . (1.1) 

In order to study the multiparticle quantum entanglement in more detail, we have to find ways how to prepare 



(synthesize) states of the form given by Eq. (1.1 



In this paper we will study in detail how N qu 



L) in various physical systems. 

ibits can be prepared in entangled states of the form (1.1). We assume 
that the qubits are encoded in internal ionic states as originally proposed in the model of quantum processor by Cirac 
. and ZoUer Q. Our paper is organized as follows: Section |l| is devoted to the description of quantum logic gates 
' and networks. Here we present multi-qubit controlled gates. We show how these gates can be expressed in terms of 
single-qubit and two-qubit gates, but we argue that for practical purposes it is more appropriate to utilize directly 
multiple-qubit gates rather than decompose them into elementary single and two-qubit gates. 



In Section III we present a logical network with the help of which symmetric states of the form (1.1) can be 



synthesized Section IV is devoted to a general problem of synthesis of a pure state of an arbitrary A^ qubit state. We 
present a simple network using which an arbitrary A'^ qubit state can be created. In Section ^ we apply this algorithm 
to a specific problem of N cold trapped ions. Following the original idea of Cirac and ZoUcr we show how the states 



of interest can be created. In the last Section VI we discuss the experimental realization of the proposed scheme 
on cold trapped ions and we also briefly address the efficiency of using multi-qubit control-NOT gates rather than a 
network of two-qubit control-NOT gates. 



II. QUANTUM LOGIC GATES AND NETWORKS 



Let us start with a brief description of those objects we will use later in the paper. We will follow the notation 
used in Ref. ij. The qubit (quantum bit) is a quantum two-level system in which logical Boolean states and 1 
are represented by a pair of normalized and mutually orthogonal quantum states labelled as |0) and |1). These two 
states form a computational basis and any other pure state of the qubit can be written as a coherent superposition 
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= a|0) + j3\\) with complex amplitudes a and /3, such that |ap + = 1. We may represent a state of a qubit 
as a point on the Bloch sphere with the parameterization a = cosi?/2 and [3 = e"''sint?/2. In quantum or atomic 
optics the qubit is often represented by a two-level atom (ion) with two selected internal levels denoted as \g) and 
|e). The quantum register of size iV is a collection of N qubits. The quantum logic gate is a quantum device which 
performs a unitary operation on selected (target) qubits conditioned by states of control qubits during a given interval 
of time. A gate acting on a single qubit is termed as a single-qubit gate, gates acting on more qubits are referred 
to as multi-qubit gates. The quantum logic network is a quantum device consisting of several quantum logic gates 
synchronized in time. 



A. Single-qubit rotation 



A single-qubit gate corresponds to a unitary operator W represented in the computational basis {|0), |1)} by the 
matrix 
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(2.1) 



A special case of a single-qubit gate is a single-qubit rotation O [see FIG. ^ (a)]. Its parameterization depends on the 
choice of coordinates on the Bloch sphere. We will define it in the matrix form in the basis {|0), |1)} as follows 



O(0,, 
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Rii j 


-i 



cos(6l/2) e*'^sin(6l/2) \ 
-e-*"*sin(6'/2) cos(6'/2) ) 



(2.2) 



where 9 refers to the rotation and <j> to the relative phase shift of the states |0) and |1) in the corresponding Hilbert 
space. 



(a) 



(b) 



(c) 



FIG. 1; A sche mati c representation of (a) a single-qubit rotation defined by the relation (2.2), (b) a two-qubit CNOT gate defined by the 
transformation (2.3) and (c) a multi-qubit (control)''-NOT gate defined by the transformation (2.4). 



B. Two-qubit and multi-qubit control-NOT gates 



A two-qubit control-NOT (CNOT) gate acts on two quantum bits denoted as the control and the target qubit, 
respectively [see FIG. |l] (b)]. If the control qubit (mi) is in the state the state of the target qubit (to2) is flipped. 
Otherwise, the gate acts trivially, i.e. as a unity operator 1. We may characterize this gate with the help of the truth 
table 



|0)mjl)m. 
|l)mJO)™, 

|l)mi |l)r?X2 



|0)™J0) 

|0)™Jl)m,, 
|l)mi |l)m2 ) 

Il)m,l0),„,, . 



(2.3) 



A multi-qubit control-NOT (CNOT) gate is defined analogically [see FIG. ^ (c)]. The only difference is the number 
of control qubits. In other words, a multi-qubit (control) '^-NOT gate acts on g -I- 1 qubits with q control qubits 
(mi, . . . , mq) and the m^+i qubit is target. If all control qubits are in the state then the state of the target qubit 
is flipped. Otherwise, the gate action is trivial. The truth table of the multi-qubit (control)*-NOT gate acting on 
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mi, . . . , TTig+i qubits reads as follows 

l^no) |0)m<,+ i 



l*no>|l)m,+ i 



l*i;e.)|0) 



14- 



yes, 



I l)iTi,+i 



l*no)|0) 



l^no)|l)m,+i , 
|^'ye.)|l)m,+ i , 
|«'ye.)|0)^,+ , . 



|*„o) ^ n '^|l>m, , 



(2.4) 



C. Multi-qubit control- i? gates 



A multi-qubit (control) "^-i? gate acts on g + 1 qubits. The mi, . . . , m^ qubits represent the contr ol p art of the gate 
while the m^+i qubit represents the target [FIG. ^. This gate performs a single-qubit rotation (2^) on the target 
qubit if all control qubits are in the state |1). Otherwise, it acts trivially. Speaking precisely, if all control qubits 
(mi , . . . , niq) are in the state 1 1) , then the operation R = R\ a R^a R2 Ri is applied (from right to left) on the 



(target) qubit. In the basis of the target qubit {|0)m^^^, 

e'^'^ sin e 



^1} we can introduce the matrices 
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-sin(6l/2) cos(6i/2) 
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(2.5) 



where i?i = 0(7r, cj)), rI = O^tt, cj)), R2 = O{0, 0) and rI = 0{e, 0). The operation 0(9, (p) is defined by the relation 



(2.2). The matrix cr denotes the NOT operation. If not all control qubits are in the state |1), then the gate performs 



on the target qubit the operation 1 = R[l R2I R2 Ri, where 1 is the unity operator. We may write the truth table 
of the multi-qubit (control) '^-i? gate as follows 



l*no)|0) 



l«'.e.)|0) 
l*ye.)|l)m,+i 



nig+l 



l*«o)|0),„,^, , 
|«'no>|l)m,+i , 



COs9\0)mq+i 

sin0|O)„,+, 



■cos 61 ID 



|l)m,+i) I 



(2.6) 



where \'^no) and |^j,es) are defined in (2.4) 
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FIG. 2: A scheme of a multi-qubit (contro^'-NOT gate acting on q + 1 qubits with q control qubits (mi, . . . ,mq). The rrto-i- ith qubit 
is the target. The operators R\,R2, R\ and R\ are defined by Eq. (2.5) and the gate is determined by the transformation ( ^.(| ). The 



gate corresponding to Eq. (2.7) is represented by the same network except the single-qubit rotations R\ and r\. 



If the preparation of a particular class of quantum states does not require the introduction of a relative phase shift 
(j) between the basis states |0) and |1), then a reduced quantum logic network is sufhcient. In particular, the operation 
R = aRlaR2 on the target qubit (m^+i) conditioned by the state of control qubits (mi, . . . ,niq) can be realized 
according to the following truth table 

.>|0)™,+, |*„o)|0) 

. l*no)|l) 



I*. 
If. 
I*. 
I*, 



9+1 ' 

)m,-|-i J 

cos6l|0),„,^i 
sin6'|0). 



■Iq + l 



- sml 
cos I 



'ID. 



(2.7) 
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The results given above for the muhi-qubit control- i? gates are compatible with the scheme proposed in Ref. [ pT[ , 
where a decomposition of multi-qubit CNOT gates into a network of two-qubit CNOT gates has been presented. 
However, this decomposition may require many elementary operations. It seems to be more appropriate for some 
practical implementations of quantum computing (for example, computing with cold trapped ions [|l6| ) to implement 
directly multi-qubit CNOT gates. 



III. QUANTUM LOGIC NETWORKS FOR THE STATE SYNTHESIS 



In this Section we present quantum logic networks for the synthesis of specific types of coherent superpositions 
of multi-qubit quantum states. Later we will use this result for construction of an algorithm for a generation of an 
arbitrary pure quantum state of N qubits. 

We will consider a quantum register of size N, i.e. N qubits. Let us denote 



Firstly, let us consider a simple network consisting of a multi-qubit control-i? gate having (N — 1) control qubits 
(ci, . . . , cat-i) and a single target qubit (t\) [see FIG. pi. Let us assume that all qubits have been initially prepared 
in the state i.e. the whole system is in the state |1) and the gate realizes the operation 

|l)^-.i?oi|l)'^^^|0)t,+i?n|l)^, (3.2) 



where i?oi and are defined by the relation (2.2). 



c 



fj \R 



FIG. 3: The network realizing the transformation given by (3.2). 



Secondly, let us consider a network with (A^ — 2) control qubits (ci, . . . ,CAr_2) and two target qubits (t\,t2) [see 
FIG. ^. The network acts on the initial state as follows (each arrow in the figure corresponds to an action of a 
gate in the sequence) 



^ i?oi|l)'^-'|0)t,+i?ii|l)^ 

^ i?oi|l)'^-'|0)tjO)t,-Fi?n|l)^-i|0), 

JV-2|n\ |n\ I D 



i?oi|l)^'-10),jO)t,+i?n|l)^V (3.3) 




FIG. 4: The network reahzing the transformation (3.3). 

Further, we design a network with (N — 3) control qubits (ci, . . . ,CAr_3) and three target qubits (ii,i2,i3) [see 
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FIG. ^. This network acts as follows 
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FIG. 5: The network realizing the transformation (3.4) 



The arrangement of quantum logic networks with more target qubits is straightforward. One has to add another 
multi-qubit CNOT gate acting on the added target qubit and then one more multi-qubit CNOT gate must be included 
at the end of the network in ord er to erase "unwanted" changes on all other terms in a superposition state [for instance 
see the 4*^* and 5*^ line in Eq. (U)]. 

As an example let us consider a network that prepares a pure symmetric (wit h re spect to permutations) entangled 
state with just one qubit in the state |0) and all others in the state |1) [see Eq. It can be shown that this state 

exhibits t he m aximum degree of entanglement between any pair of N qubits ]14|. The network for the synthesis of 
the state (1.1) from the initial state |1)^ is shown in FIG. |[ where the rotations Uj are defined as follows 



1 / 



3 



,iV- 1 



(3.5) 



The action of the network in FIG. can be described as follows 
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-1= V |1)^-1|0), + ^|l)--l|0)Ar = ^ V |1)^-^|0), 
.7 = 1 .7 = 1 



(3.6) 



where |1)— denotes the state with all qubits in the state |1) and |1)— ^^|0)j rep rese nts the state of the register with 
{N — 1) qubits in |1) and the jth qubit in the state |0) [see the notation in Eq. (3.1)]. 

A very simple example is the synthesis of the GHZ state, i.e. a coherent superposition with all qubits to be in the 
state JO) or |1) with the same probability, i.e. \'E.)ghz ~ (|0 )^ + |l)^)/-\/2. The corresponding network is shown in 
FIG. 0. The single-qubit rotation R = 0{tt/2,tt) defined in (2.2) is applied on the initial state |0)— and prepares the 
superposition (|0)- + |0)-^i|l)i)/y2. Applying sequentially aU CNOT gates one prepares the GHZ state \^)qhz- 
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FIG. 6: The network for the synthesis of the symmetric entangled state (1.1) on A'^ qubits. The rotations Uj are given by Eq. (3.5). The 
A' qubits are assumed to be initially in the state |1)'^. 
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FIG. 7: The network for the synthesis of the generalization of the GHZ state. 
R = 0{-k/2,tt). The initial state is |0)^. 



The single-qubit rotation _R is given by Eq. (2.2) for 



IV. SYNTHESIS OF AN ARBITRARY PURE QUANTUM STATE 

Coherent manipulation with states of quantum registers and, in particular, the synthesis of an arbitrary pure 
quantum state is of the central importance for quantum computing. One of the important tasks is the preparation of 
multi-qubit entangled states. 

Based on the discussion presented above we can propose an array of quantum logic networks that prepare an 
arbitrary state from the register initially prepared in the state |0)^, i.e. 

2"-l 11. ..1 

io)~ mN))^ E ^^■1^^)= E ^-1^)' (4.1) 

j=o £!;=00...0 

where a; is a binary representation of the number 2^ . The proposed scheme can be generalized on the quantum register 
of an arbitrary size, but for simplicity we will firstly consider the case of three qubits. 
A general state of three qubits is given as 

|V;(3)) = ao|000) + e''^iai|001) + e^'^=a2|010) + e''^^a3|100) 

+ e*'^''a4|011)+e*'^5a5|101) + e*'^«a6|110)+e*'^^a7|lll), (4.2) 

where ao; ■ • ■ > ci7 are real numbers satisfying the normalization condition 
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E«- = l, (4.3) 

3=0 

and ifii, . . . ,ip^ are relative phase factors. The global phase is chosen such that = 0. 

In what follows we will prese nt th e procedure for the synthesis of the state (4^). Let us use the abbreviated form 
of the matrix R defined in Eq. (2.5) which we denote as 
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FIG. 8: An array of networks for the synthesis of an arbitrary pure quantum state (4.S) on three qubits. The initial state is |000) and 
the rotations Vj are given by Eq. (4.4). 
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FIG. 9; A compact form of the array of the networks shown in FIG.|8| 



where = cosBj and bj — sin0j. The initial state is |000). The network presented in FIG. || (a) prepares out of the 
state 1 000) the superposition 



ao|000)-e-*2'^'o5o|iii) 
Applying the network in FIG. § (b), a new term 

_g»2(0i-0o)5j^^^|OOl) 



(4.5) 



(4.6) 



is added to the superposition (4.5) while the amplitude of the component |000) is not affected at all. The application 
of the network given by FIG. |^ (c) adds another new term 



*2(02-0o)5pQ^52|OlO) 



(4.7) 



and does not influence the amplitudes of two foregoing terms |000) and |001). Repeating this procedure, the network 
in FIG. H (d) adds a new term 



_g»2(03-0o)5p^j^Q^Q253|lOO). 

Analogously, the network shown in FIG. |^ (e) adds a new term 

-6*2(04-00)5^^^^203641011). 
while the networks (f), (g) shown in FIG. || (f) and (g), add new terms 

_g^2(05-0o)5g(j^a2a3a465|lOl), 



(4.E 



(4.9) 



(4.10) 



- e*2('^«-'^'')6oaia2a3a4a566|110), 



(4.11) 



8 



j 






state 





Co 


(default) 


000 


1 




2(61 - 6q) +n 


001 


2 




2(02 - 5!>o) + TT 


010 


3 




2(<A3 - (An) + TT 


100 


4 




2(04 - (?!)o) +7r 


Oil 


5 




2(05 - 0o) +vr 


101 


6 


6oaia2a3a4a5&6 


2(06 - 0o) +7r 


110 


7 


60111203040506 


-200 + TT 


111 



TABLE I: The network in FIG. ^ generates the state (4.2) from empty register |000). The network is characterized by the coefficients 
a,- , (A,- , where bj = The state ( ^ is determined by the coefficients aj,ipj. The table relates these two set of numbers. The 

inverse relations are given by the equations (4.13) and (4.14). 



respectively. The last network shown in FIG || (g) also determines the amplitude of the last term 

- e^2W7-<^o)^gQ^a2a3a4a5a6|lll) . (4.12) 



Comparing the output from the networks shown in FIG. determined by the relations (4.5)-(4.12), with the expression 
(U), we get the final results in TABLE |. 



The coherent super pos ition ( |4.2| ) is completely determined by 15 parameters (ao, . . . , 0:7; (pi, . . . , (py). The nor- 
malization condition ( fl.3| ) reduces this number to 14. The n etw orks in FIG. || are determined by 14 parameters 
(bo, ■ ■ ■ ,be; (f)Q, . . . , (po). Thus, the mapping between the state ( |4.2| ) and the networks is clearly defined. From given 
values of aj and i^j one can calculate hj and (pj according to the expressions 



and 



i(7r-937), = i(^j- - (^7) , j = l,...,6 (4.13) 



bo^ Jl-al, ^ , =, j = l,...,6, (4.14) 




which determ ine the single-qubit rotations (4.4) 



The state (4.2) contains terms corresponding to all possible permutations of three qubits. However, a reduced 
superposition with some terms missing might be desired. For this purpose, we can skip networks responsible for the 
synthesis of these terms or the corresponding parameter hj can be set to zero. For instance, in the case when the term 
|000) does not appear in a final desired quantum state, we begin with the initial state |111) and skip the network in 
FIG. H (a). If we do not wish, for a change, to generate the term |111), one may set the parameter ag to zero and the 
phase factor can be chosen arbitrarily (see the table above). 

The scheme can be analogically extended to an arbitrary number of qubits. In what follows we will briefly discuss 
the extension on four qubits. These can be prepared, in general, in the coherent superposition consisting of 16 terms, 
i.e. 10000), 10100), |0010), ... , |1111). 

The network in FIG. |l^ (a) prepares the superposition of the terms |0000) and |1111) with correspond- 
ing complex amplitudes depending on the choice of the single-qubit rotation Application of the net- 
work in FIG. M (b) running through all possible permutations of four qubits, i.e. (ci, ii, t2, ^3) = 
{(1,2,3,4);(2, 1,3,4);(3, 1,2,4);(4, 1,2,3)}, adds to the superposition new terms 1 1000), |0100), |0010) ,10001) with 
corresponding amplitudes determined by i?2- Further, we apply the network of the type in FIG. M (c) run- 
ning through the permutations (ci, C2, <i, ^2) = {(3, 4, 1, 2); (2, 4, 1, 3); (2, 3, 1, 4); (1, 4, 2, 3); (1, 3, 2, 4); (1, 2, 3, 4)} and 
the terms |0011), |0101), |0110), |1001), |1010), |1100) (with correspondingamplitudes given by R^) wiU be in- 
cluded to the state under construction. Finally, the network in FIG. |o| (d) running through (ci, C2, C3, ti) = 
{(2, 3, 4, 1); (3, 4, 1, 2); (4, 1, 2, 3); (1, 2, 3, 4)} generates new terms |0111), |1011, |1101), |1110). 

The extension to N qubits is analogical. The state synthesis is started from the initial state |0)^. Firstly, one 
uses the network for the preparation of superpositions of |0)^ and |1)^ with determined amplitudes. Secondly, the 
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FIG. 10: An array of networks for the synthesis of an arbitrary pure quantum state (1.1) of four qubits as discussed in Section |lV[ ) . 



networks with one control qubit (ci) and iV — 1 target qubits (ii, . . . ,iAr-i) running through all permutations are 
applied. Then, we employ the networks with two control qubits (ci,C2) and N — 2 target qubits (^i, . . . , iAr_2). 
Further, the networks with more control qubits (3,4, . . . ,iV — 1). These procedures are repeated until we achieve 
A'^ — 1 control qubits (and one target qubit). The synthesis stops and a desired final state is prepared. 



V. REALIZATION ON COLD TRAPPED IONS 

In previous sections we have proposed a scheme for the synthesis of an arbitrary pure quantum state of a system 
of qubits. The implementation of the multi-qubit CNOT gate has played the central role in our scheme. It is well 
known how to decompose multi-qubit gates into a network of single-qubit and two-qubit CNOT gates However, 
it seems that a direct implementation of multi-qubit CNOT gates in specific quantum systems is more straightforward 
and requires less elementary operations (for example laser pulses) than its decomposition. 

We demonstrate this idea on a system of cold trapped ions. We will briefly describe the system under consideration 
and show how multi-qubit gates can be implemented. 

The quantum system considered here is a model of a string of N atomic ions confined in the linear Paul trap 
proposed by Cirac and ZoUer in 1995 ||l^. First experiments on a single ion and two ions were realized by the NIST 
group in Boulder | [T^ . Experiments with more ions were done, for example, by the group in Innsbruck jigf . 

The confinement of a system of trapped ions along the a;, y and z axis can be described by an anisotropic harmonic 
pseudopontential of frequencies lo^ <C Wr, where for the usual choice of trapping radio- frequency (rf) voltage we get 
The ions are firstly Doppler cooled and then undergo the sideband cooling. Laser cooling minimize 
their motional energy and the ions oscillate around their equilibrium positions. In this case we can describe their 
motion in terms of normal modes. We will consider only the lowest, center-of-mass (COM), vibrational collective 
mode of the ions along the z axis, when all the ions oscillate back and forth as if they were a rigid body. The sideband 
cooling leaves the ions in the quantum ground motional state, therefore we have to assume the Lamb-Dicke limit, 
i.e. the photon recoil frequency (corresponding to the laser cooling transition) is much smaller than the frequency of 
the considered COM mode. The ions in the trap represent qubits with two distinct internal atomic states denoted 
as \g) and |e) with corresponding energy levels Eg and -Eg, respectively. We will consider individual- ion-addressing 
with a laser beam of the frequency ul represented by a classical traveling wave. Then, in the interaction picture, 
in the rotating wave approximation plus the weak coupling regime and in the Lamb-Dicke limit we can write the 
Hamiltonian corresponding to the interaction between the jth trapped ion {j = 1, . . . , N) and the laser beam tuned 
on the carrier [uj^ = ^o) 



and on the first red sideband [lol — i^o ■ 



irj 



|e)j(ff|a+ l5)j(e|a^ 



(5.1) 



(5.2) 
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where flj = \i}j\e~'^'^ is the laser couphng constant, </> is the laser phase, rj is the Lamb-Dicke parameter, a and are 
the annihilation and creation operator of the quantized COM mode with the frequency cu^, where a^a\n) — n\n) and 
ujQ — {Ee — Eg)/'h is the atomic transition frequency. 

Further, we can write the unitary evolution operators via which the action of quantum gates is realized. Firstly, 
let us consider the evolution operator corresponding to a fcTr-pulse on the carrier {t — kTT/\D,j\) applied on the jth ion 
with the arbitrary initial choice of the laser phase such that 

kn 



4(0) = exp 



191 e 



J9 



Under the action of this unitary operator the two internal states of the j-th ion are changed as follows 

\g), cos(fc7r/2)l5>j--e"**sin(fc7r/2)|e),, 
\e)j ~* cos(fc7r/2)|e)j +e*"*sin(/c7r/2)|5)j . 



(5.3) 



(5.4) 



Secondly, we have the evolution operator for a fcTr-pulse on the first red sideband {t = hnV N /\D,j\r]) on the jth ion 
choosing the laser phase such that 



i3^«(0)=exp 



(5.5) 



which implies the transformation 

\9)j\l) 



l5>,|0) 



,|0) 



cos(fc7r/2)|g)j|l) - ie"^"^ sin(fc7r/2)|e)j|0) , 

cos(fc7r/2)|e)j|0) - ie*"^ sin(A:7r/2)|5)j|l) , (5.6) 

where q — I,II and |ej) denotes the upper internal level, whereas \ejj) refers to an auxiliary internal level \aux). In 
the original proposal llq] the values of the parameter q = I,II refer to the situation where the transition excited by 
the laser depends on the laser polarization. 

The operators (5.3) and ( |5.5| ) provide us with the possibility to introduce the implementation of the single-qubit 
rotation and multi-qubit CNOT gate on selected ions (representing qubits). It is obvious from the transformation (5.4) 
t hat the evolution operator ( ^.3D corresponds to the single-qubit rotation 0(fc7r, (p) on the jth ion [sec the definition 
(2.2)]. The two-qubit CNOT gate (the mith ion is the control and the m2th ion is the target) is realized by the 
evolution operator (from right to left) 

iV2(,)e^/6^//e^/il/2(0), 

which corresponds to a sequence of pulses as described above. This transformation acts on two ions as 

l5)mi|5)m2|0> ' \9)ra,\9)mJO) , 



(5.7) 



\9)r 



^l|e)m2|0) 
Ill5)m2|0) 



\9)r 



JO), 
JO), 



|e)™Je)„jO) |e)™j5)™jO). (5.8) 

The ions are assumed to be cooled to the ground vibrational state |0) before the operation. We have used the notation 
13 = B{0) in the relation (5.7). The two-qubit CNOT gate can be extended to the multi-qubit (control) ^-NOT gate 
acting on q + 1 ions (toi, . . . , ions represent the control, while the mg+ith ion is the target) and can be realized 
by the following evolution operator (from right to left) 



corresponding to the transformation 

I*«o)|g)™,+J0) 



Ie)™,+J0) 



yes)\9)ra,+ M 
yea) \^)mq-^i |0) 



>l5>m,+ JO), 

>|e)m,+i|0), 
Ne)™,+J0), 

yes/ |ff)m<,-|-i |0) . 



I*,: 



(5.9) 



yes/ 



I* 



yes/ 



n 



(5.10) 
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It is obvious from Eqs. (5.8) and (5.10) that the ions must be kept in the ground motional state. This arrangement 
eliminates heating processes which lead to decoherence. However, it is still the experimental challenge to cool to the 
ground state |n = 0) more than two ions. 

VI. DISCUSSION AND CONCLUSIONS 

In this paper we have shown how multiparticle entangled states can be constructed with the help of multi-qubit 
quantum gates. We have shown how to implement these gates on the system of cold trapped ions. This allows us to 



"realize" any multi-qubit control- i? gate and also any logic network proposed in Sec. Ill and IV. 

To understand the feasibility of this algorithm we present some estimations considering the application of the 
introduced gates and networks on cold trapped ions. 

The main aim of further discussion is to illustrate a range of relevant physical parameters for implementation of 
proposed scheme. Obviously, specific experimental setups have to be considered separately. We present just rough 
estimates of minimal times required for realization of desired gate operations. 

Let us consider Calcium ions ''"Ca^ with the "ground" (computational) state \g — S1/2) and the "excited" (com- 
putational) state |e = 0^/2)- The lifetime of the ion on the metastable -D5/2 level is 1.045 s. 

We will assume N ions loaded and confined in the trap. The ions will be individually addressed with a laser beam 
(A = 729nm) supposing the Gaussian intensity profile I / lo = exp(— 2/32/i«g), where p denotes the radial distance and 
2wq = 10/im is the beam waist. Further, let the angle between the laser beam and the z axis he d — 60°. Then, the 
recoil frequency of the Calcium ion is Jr = 2.33kHz, where Jr = ER/h, Er = h'^k'^/2m, k = 27r/A and h = 2T:h. The 
axial trapping frequency is uJz/^tt = llOkHz. We can also calculate the Lamb-Dicke parameter 77 = ^jERjhljzi i-e. 
r] = 0.15. The minimum spacing between two neighboring ions is determined by the approximate formula pQ, pl| : 

2.018 / V'' 

where q is the ion charge, m is the ion mass and eo is the permitivity of vacuum. 

The multi-qubit CNOT gate on the ions is realized by the evolution operator (5.9) . W e will consider three types 
of elementary operations: (1) 7r/2-pulse on the carrier [A) defined by the relation (|5.3|), (2) 7r-pulse (B^) and (3) 
27r-pulse (S^) on the first red sideband ( |5.5| ). Each elementary operation takes a certain time to be implemented on 
the system of cold trapped ions. Steane et al. addressed in detail the speed of ion trap information processors in [p2[ . 

Firstly, the single-qubit rotation [A) can be made much faster than two-qubit operations (S^,^^), because the 
Lamb-Dicke parameter 77 can be set zero (i.e., the laser beam is perpendicular to the z axis). Thus, can be made 
large without restrictions on the weak coupling regime characterized by the condition <C oJz- We will assume 
|il|/27r = 50kHz and estimate the time required for the single-qubit rotation as T4 = 7r/2|ri| = 5/iS. 

Secondly, by definition for the operations and , the Lamb-Dicke parameter must be non zero [see Eq. ( ^.2[) ]. 
This means that some unwanted off-resonant transitions will be present, which may significantly affect times required 
for the operations B^'^ 

In Ref. p2| it has been shown that the minimal time Tg for the realization of the operation B^ is proportional to 
the geometric mean of the recoil and trapping frequency, i.e. 



1 2^/2e IEruj^ 



Tb ^ \ h 2n 



(6.2) 



where the imprecision e = \/l — F is defined via the fidelity F of the process. The time for the operation B^ is then 
2Tb. 

Once the gate times are estimated, we can determine the minimal total time T required for the experimental 



preparation of the state (1.1) on Calcium ions. The total time T is the sum of times of all operations A^ B^,B^, which 



appe ar in the implementation of the network in FIG. |^. The total number of all operations, when preparing the state 



(1.1) on N ions, is 2N'^ + AN — 10. The explicit expression for the total time reads 

T ^ N{A)Ta + N{B^)Tb + N{B'^)2Tb. (6.3) 

In what follows we will consider several situations with the number of trapped ions varying from 2 to 20. In a 
given ion trap for different values of ions we obtain different minimal spacings Azmin [see Eq. (^^)]. The minimal 
spacing between ions has to be larger than the half-width of the Gaussian profile of the addressing laser beam. In 
the Innsbruck experiment ]2^ the width of the Gaussian profile is proportional to 10/zm. Even for 20 ions with 
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Azmin = 7.31^m [see Eq.(6.1)] and the given width of the Gaussian profile, the ratio between the hght intensity of the 
laser addressing a given ion to the intensity of the same beam on the neighboring ion is as small as 1.4%. Therefore, 
individual ions can be addr essed rather efficiently. 

As follows from Eq. (x2) the minimal time for the gate operation depends on the required fidelity of the process. 
In our case we consider two values of the fidelity, namely F — 99% and F = 75%. Given these values we can estimate 
relevant physical parameters. 

In TABLE ^ we present results of our estimations. From here we can conclude that for a given lifetime of Calcium 
ions (1.045 s) one can perform in our scheme a coherent manipulation with up to 20 ions with the fidelity 99%. It 
seems to be a very optimistic estimation, however we did not optimize the network itself. 



N 




Tb 


/is] 


N{A) 


iV(S') 




T [ms] 


F = 99% 


F = 75% 


F = 99% 


F = 75% 


2 


24.4 


312 


62.4 


3 


2 


1 


1.26 


0.265 


3 


20.8 


382 


76.4 


9 


8 


3 


5.39 


1.11 


4 


18.0 


441 


88.3 


15 


18 


5 


12.4 


2.55 


5 


15.9 


493 


98.7 


21 


32 


7 


22.8 


4.65 


6 


14.3 


540 


108 


27 


50 


9 


36.9 


7.48 


7 


13.1 


584 


117 


33 


72 


11 


55.1 


11.2 


8 


12.2 


624 


125 


39 


98 


13 


77.6 


15.7 


9 


11.4 


662 


132 


45 


128 


15 


105 


21.1 


10 


10.8 


698 


140 


51 


162 


17 


137 


27.7 


15 


8.59 


855 


171 


81 


392 


27 


382 


76.7 


20 


7.31 


987 


197 


111 


722 


37 


786 


157 



TABLE II: N is the number of Calcium ions in the trap, Azmin is the minimal distance between two neighboring ions ( |6.l| ), Tg is the 
minimal time for the realization of the operation [in Eq. (p^)] for two different values of the fidelity (_F = 99%, F = 75%). N{jV) is the 
total number of the operatio ns A in the network in FIG. ^ N{B'^) and N{B^) a re to tal numbers of the operations B^ and B^, respectively. 
T is the total minimal time ( ^.3| ) for the experiment al p reparation of the state (1.1) on A'^ ions via the network in FIG. ^. T4 = 5/is is the 
time for the realization of the operation A [in Eq. (p.9[)]. 



We have chosen the cold trapped ions as an example for the situation when the direct implementation of the 
multi-qubit CNOT gate (using elementary operations, i.e. in this case laser pulses) is much less demanding than the 
decomposition of multi-qubit CNOT gates into the network of two-qubit CNOT gates. For instance, let us consider 
the multi-qubit CNOT gate on six qubits. Using results of Ref. [0 we can decompose this multi-qubit CNOT gate 
into the network composed of 12 two-qubit CNOT gates. In additi on, this network had to be extended by three 
additional auxiliary qubits. The multi-bit CNOT gate on N ions (5.9) is r ealized by 2N + 1 laser pulses. Each 
two-qubit CNOT gate on two ions is then realized using five laser pulses ( |5.7| ). It means that all together 60 pulses 
have to be used for 12 two-qubit CNOT gates. However, the direct implementation of the multi-qubit CNOT gate 
on six ions requires only 13 laser pulses. This difference becomes even more significant with the increasing number of 
the ions. Obviously, smaller the number of pulses easier the scheme can be implemented. 
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